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ABSTRACT 

We examine the expected radio confusion on the thermal Sunyaev-Zel'dovich 
(SZ) effect of galaxy clusters at 28.5 GHz and 90 GHz based on the cluster radio 
luminosity function (CRLF) at 1.4 GHz. With the observationally determined 
spectral index (a, where S v oc v~ a ) distribution, instead of a single average 
index (a), we convert the cluster radio luminosity function at 1.4 GHz to the 
high frequency ones and estimate the total radio flux in a cluster. At 28.5 GHz, 
radio confusion is up to 10 ~ 100% for small clusters (M ~2x 1O 14 M ) below 
redshift z=l, with more severe confusion for smaller and lower redshift clusters. 
By contrast, at 90 GHz the confusion is less than 10% for small clusters even at 
low redshifts. 

Subject headings: cosmology: theory — galaxies: clusters: general — radio con- 
tinuum: general 

1. Introduction 

The cosmic microwave background radiation (CMBR) provides a light source for ex- 
ploration of the structure formation history. Before the CMBR photons reach the present 
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epoch, they interact with the hot diffuse ionized gas residing in the galaxy clusters over a 
wide range of redshifts. This process results in subtle changes to the CMBR spectrum due 
to the inverse Compton scattering, the so-called Sunyaev-Zel'dovich (SZ) effect (Sunyaev & 
Zel'dovich, 1972; see Birkinshaw 1999 and Carlstrom et al. 2000 for reviews). In the past 
two decades the detection of SZ effects has been reported at an increasing rate (Grego et 
al. 2000; Carlstrom et al. 2000; Joy et al. 2001; Grego et al. 2001) and this effect can 
therefore provide a realistic powerful probe for cosmological and cluster studies. Since the 
SZ flux is not affected by the distance but depends only on the intrinsic thermal energy of 
hot gases contained in the system, the SZ observation, in conjunction with observations of 
other wavelengths, such as X-ray emitted from intracluster gas, and optical weak lensing 
measurements, enables us to determine the cosmological parameters by studying the high 
redshift galaxy clusters. The SZ blank sky surveys have been suggested (Bartlett & Silk 
1994; Barbosa et al. 1996). Several works follow (Bartlett 2000; Holder et al. 2000; da Silva 
et al. 2000; Kneissl et al. 2001; Fan & Chiueh 2001; Xue & Wu 2001), and have predicted 
how the SZ cluster counts can be related to cosmological parameters. In addition, several 
next generation interferometric arrays are either under construction or have been proposed, 
e.g., AMI (Kneissl et al. 2001), SZA (Carlstrom et al. 2001), and AMiBA (Lo et al. 2001). 

Radio point sources or radio galaxies are usually found in the galaxy clusters. Therefore, 
the accuracy of the millimeter to centimeter wavelength observations of the SZ effect can 
be contaminated by the emission from radio sources (Loeb & Refregier 1997). Below 218 
GHz, at which the thermal SZ effect becomes null, one measures the SZ intensity or flux 
decrement. The radio sources, however, produce excess emission so that one can underesti- 
mate the measured SZ decrement, yielding systematic errors in the estimate of cosmological 
parameters. This confusion has been estimated by Cooray et al. (1998, hereafter CGHJC), 
who observed toward Abell 2218 at 28.5 GHz and concluded that the correction to the Hub- 
ble constant is less than 6%. In addition to point sources, diffuse cluster radio sources, i.e. 
halos and relics, (e.g., Giovannini et al. 1993; Rottgering et al. 1997) also lead to confusion 
on the cluster SZ. 

With the requirement of accurate SZ maps, the radio source subtraction becomes essen- 
tial for the low brightness (order of 100 /zJy and mJy) SZ measurements. For this reason, 
recent SZ observations are performed in combination with the measurements of radio sources 
(CGHJC; Komatsu et al. 1999; Pointecouteau et al. 2001). While one may remove strong 
sources as they can be detected above the flux limit, the faint sources below the flux threshold 
blend with the SZ signal, and result in inaccuracy. However, there have not been systematic 
surveys on the radio sources in high redshift clusters nor systematic surveys at frequen- 
cies higher than 30 GHz. In fact, the high-frequency (e.g. 90 GHz) radio survey at the 
mJy flux limit is a formidable task, as the present telescopes are not sufficiently efficient 
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in conducting an unbiased survey with a wide field coverage and high sensitivity. In view 
of the lack of reliable information for estimating the radio contamination on the SZ clus- 
ter surveys at high frequencies, such as AMiBA (90 GHz), we are motivated to derive an 
empirical high-frequency cluster radio luminosity function (CRLF) by synthesizing available 
low-frequency data (Ledlow & Owen 1996, hereafter L096; Slee et al. 1996, hereafter SRA96, 
and CGHJC). Unfortunately, the spectral evolution of radio sources is rather uncertain. Due 
to the uncertainties of redshift evolution, it is reasonable, as a first attempt, to adopt a no- 
evolution hypothesis to evaluate the radio flux in high redshift clusters. In this paper we aim 
to calculate the expected confusion of high-frequency radio point sources on the cluster SZ 
effect over substantial redshift and cluster mass ranges. In section 2, we briefly review the 
formulation of the SZ flux. Section 3 predicts the high frequency cluster radio luminosity 
functions (CRLFs) based on low frequency surveys. The expected high frequency radio flux 
are calculated in section 4. Finally, the conclusion is given in section 5. Throughout this 
paper, we use H = 65 km s -1 Mpc -1 , = 0.35, and fl\ = 0.65. 



One may evaluate the total SZ flux by the procedure as follows (see also Barbosa et al. 
1996, Xue & Wu 2001). For the CMB radiation passing through the intracluster hot gas, 
the change in CMB intensity by thermal SZ effect is 



where m e is the electron rest mass, n e is the electron number density, T e is the electron 
temperature, and <tt is the Thomson cross section. j„(x) represents the spectral dependence 



in which x = hpu/kTo and To=2.728 K is the CMB temperature. For v = 90 GHz, x ~ 1.58. 
The total SZ flux of a cluster at redshift z is the integral of the intensity over the solid angle 
subtended by the cluster 



2. 



SZ Flux 




(1) 






(3) 



where D A is the angular diameter distance to the cluster. If we assume the gas is isothermal, 
then the total SZ flux depends only on the electron number, which can be replaced by the 
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gas mass of the cluster. Introduce the gas fraction / b = M gas /M vir , where M vir is the virial 
mass of the cluster, then the equation above becomes 

D z A (z) \m c c 2 J \/j, c m p J 

where /x e = 2/(1 + X), and X=0.768 is the hydrogen mass fraction in the primordial abun- 
dances of hydrogen and helium. 

Define i? v ; r as the radius of the cluster within which the mean density is A c times the 
critical density of the universe p c at a certain redshift, i.e., M vir = 4nRl ir p c A c /3, then the 
relation between temperature and mass is (Bryan & Norman 1998) 

kTe = L39/t {i^k ) 2/3 ^ AcEiz) ^ 1/3 keV ' (5) 

where E(z) 2 = f2 M (l + z) 3 + f2 K (l + z) 2 + f2A, = 1 — — ^a, and the redshift-dependent 
Hubble constant H(z) = 100hE(z) km s _1 Mpc -1 . / T is a normalization factor obtained 
by numerical simulations and about 1. A c is the density contrast of the virialized spherical 
halo to the critical density of the universe at that redshift, which can be fitted as (Eke et al. 
1996) A c = 18tt 2 + 82[Q(z) - 1] - 39[Q(z) - l] 2 for fi K = 0, where Q{z) = Q M (1 + zf/E 2 (z). 

The gas fraction f h varies with the cluster temperature by the empirical formula (Mohr 
et al. 1999) 

/ hT \ °- 34±0 - 22 

/ bi50 = (0.207 ± 0.011) t^—) , (6) 

where the subscript 50 denotes h =0.5. We may scale /b by /b = fb,5oh^' 5 (White et al. 
1993), and Eq.(4) and (6) will be used to compare with the radio flux calculated in the next 
section, so as to determine the radio confusion on the SZ flux. The reason for incorporating 
the temperature dependence of f h is as follows. It is likely that preheating by supernovae 
and/or AGNs provides nongravitational heating to the intracluster gas, which makes the 
distribution of the hot gas extend to large radii (Wu et al. 1998, 2000). The effect is more 
significant in poor clusters than in rich ones because the former have a shallow gravitational 
potential, as evidenced by the steepening of the L x -T relation (David et al. 1993; Wu et al. 
1999) as well as the excess entropy in galaxy clusters (Ponman et al 1999). 



3. CRLF and Radio Flux 

To estimate the total flux of the radio galaxies in clusters, we adopt the CRLF given by 
L096 (Paper VI). They did a series of surveys on a large sample of radio galaxies in Abell 
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clusters using VLA, e.g., in Paper IV (Ledlow & Owen 1995a, hereafter L095a), where the 
radio sample and cluster properties were presented; in Paper V (Ledlow & Owen 1995b), they 
did optical observation and investigated the optical properties toward their radio sample. 

In L096, the CRLF is for a complete sample of Abell clusters with z < 0.09 and can be 
fitted by a continuous broken power law 

1 t 1 ^radioQogPl^) , r , D 

log/i. 4 = log = a + 61ogPi. 4 , (7) 

jV op t 

in which N opt represents the number of galaxies over all magnitudes brighter than —20.5 
within 0.3 Abell radius of the cluster center and A^ radio (logPi. 4 ) represents the number of 
radio galaxies which emit the radio power P 1A at 1.4 GHz. The coefficients b and a are fitted 
with the central values as -0.15 and 1.77 for P 1A < 10 24 8 WHz" 1 , and -1.43 and 33.67 for 
Pia > 10 24 - 8 WHz-\ respectively. 

We evaluate the predicted cluster radio luminosity functions at other frequencies, f v , 
based on the 1.4 GHz survey by convoluting f\ A with the spectral index distribution n(a) 
as 

log/„ = J J (log fiA)n(a)S log P 1A - log P v - a log ( — ^ H J dadlogP 1A , (8) 

where f\ A is the luminosity function given by Eq.(7). The spectral index distributions n(a) 
that we adopt are fitted by two Gaussians and normalized as 

f o.21e- a72 ( Q - L14 ) 2 + 1.06e- 110 ( Q - a98 ) 2 ; SRA96 sample 
n[a) ~ { O.^e-^- 1 - 02 ^ + 0.44e- 3 - 41 ^ - 57 ) 2 ; C GH J C sample ' 

and are shown in Figure 1. The circles are from SRA96 and the crosses are from CGHJC. 
The two distributions are normalized to their total number of sources for comparison. The 
SRA96 sample contains 254 radio sources and its n(a) appears more regular than that 
of CGHJC sample, which contains 53 sources (a available) only in very massive clusters 
(M > 1O 15 M ). Moreover, the spectral index in SRA96 sample is determined at lower 
frequencies (1.5 GHz < v < 4.9 GHz), whereas CGHJC sample at higher frequencies (1.4 
GHz < v < 28.5 GHz). Despite these differences, the two n(a) do not differ much, at least in 
the main body of the distribution. We show in Figure 2 the predicted CRLFs using SRA96's 
n(a) at v = 28.5(1 + z) GHz, and 90(1 + z) GHz for z=0.25, respectively. Note that the 
cosmological expansion has been considered, where the frequency of the predicted function 
has a redshift dependence v = uq(1 + z), where uq is observed frequency. For simplicity in 
calculating the total flux, we fit the luminosity-weighted CRLFs by a broken power law, 

W(P f) = I CiM +7iMlogP„; logP < logP b 
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where log Pb is the break point of a continuous luminosity-weighted CRLF, and the fit coeffi- 
cients are list in the Table 1. We see that both the break points of SRA96 and CGHJC vary 
slightly with frequency, and that the peak power of radio sources decreases with increasing 
frequency. But the fit coefficients from the two groups are not in good agreement, as shown 
below. This is due to that the sample of CGHJC has flatter spectra, and the resulting power 
is higher. 

Assuming neither evolution effect on the spectral indices nor on the intrinsic bright- 
ness of radio sources, we obtain the best fit for the frequency dependence of these CRLFs 
coefficients, 




and 



CM 

C 2 {y) 

72 W) 



1.68 + 0.015i/ -2.1 x 10~ 5 z/ 2 
0.85 - O.OOIOz/ + 1.91 x 10~ 6 z/ 2 
23.84 + 9.92e-°- 015iy 
-OA4e~ om7u - 2.9 x 10~ 4 u 



2.68 - 0.95e-°- 015,/ 
0.85e~ om27u + 1.49 x 10 3 ^ 
30.42 + 3.29e-° m5u 
-0A3e-° m22u - 6.14 x 10~ 4 u 



SRA96's n(a) 



;n) 



CGHJC's n(a). 



(12) 



Figure 3 shows how the coefficients depend on frequencies and their best-fit curves. One may 
construct the CRLF of any redshift by interpolating the emitted frequency u, as u (l + z), 
where v$ is the observed frequency. 

Taking an integral over the radio power range from 10 17 to 10 26 WHz -1 and multiplying 
it by the total number of cluster galaxies, we have the total radio power emitted from the 
radio galaxies within a single cluster 



ptot . 
1 v 1 v vir 



(13) 



where N vir is the number of galaxies in a cluster of mass M v - ir . We estimate the total number 
of radio sources in a cluster by the mass-to-number relation for a given cluster mass M v \ r 
(Carlberg et al. 1996). They use the CNOC Cluster Survey to derive the M vir /N ga \ ratios 
using the sample with r-band absolute-magnitude limits from M r K = —18.0 to —20.0. To 
be consistent with A^pt in Eq.(7), which has the optical limiting magnitude at —20.5, we fit 
the mass-to-number relation derived by Carlberg et al. (1996), and extrapolate the limiting 
magnitude to M r K = —20.5 to obtain M vil /N gal ~ 1.5 x 1O 13 /i _1 M . For a given cluster mass 
M vir , we take this A^i as the number of galaxies A" v i r in Eq.(13). 
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To be compared with SZ flux, the total radio flux observed at frequency v at redshift 
z can be evaluated from the total power at frequency v as 



where D-^(z) is the luminosity distance to redshift z, D-^(z) = (1 + z) 2 Dx{z), and radio 
sources are assumed to have no evolution. We combine Eq.(4) and Eq.(14) to determine the 
confusion. 



In Figure 4, we plot the total radio flux v.s. the virial mass of a galaxy cluster (or 
group) at redshift ,2=0.25, 0.5, and 1.0. The left panels are for the observed frequency 28.5 
GHz (operated in SZA) and the right for 90 GHz (operated in AMiBA). The la error bars 
result primarily from the variance of the predicted luminosity-weighted function P u f u , and 
are much larger than the observational errors of f\ A and 71(a), which are neglected. These 
large error bars are due to the broad scattering of the distribution of the spectral index a. 
As a consistency check, we compared our estimate flux at 28.5 GHz with the CGHJC sources 
in Figure 4, and at 1.4 GHz with the L095a sources in Figure 5. The dash lines are our 
predicted mean flux, and the data points are the radio flux provided by L095a. We have 
investigated the redshift and the X-ray temperature of their selected clusters (c.f. Ebeling et 
al. 1996, 1998) and converted the cluster temperature to the virial mass by Eq.(5) to obtain 
the coordinates of the data points in Figure 4 and 5. It shows that the total radio flux at 
28.5 GHz of most but few clusters in their sample are consistent with our predictions within 
lei. The same agreement is also found at 1.4 GHz. Since the measured clusters of CGHJC 
are massive, their sample may have some selection bias. It may result in asymmetry in their 
spectral index distribution 71(a), leading to the predicted radio flux higher than those from 
SRA96, by almost a constant factor ~ 2.5 — 3 throughout all frequencies of interest. It is 
instructive to examine the averaged spectral index (a) derived from each individual n(a) 
given by SRA96 and by CGHJC. It is found that 



vo,obs 



P* ot (l + z) 
AirDl(z) 



(14) 



4. 



Result 




-1.15; forSRA96 
-0.9: for CGHJC 



(15) 



The (a) value for CGHJC is somewhat higher than that for SRA96, but is lower than the 
(a) value, —0.77, given by CGHJC based solely on their 53 sources. The flatter spectrum 
signifies the somewhat different radio sources in massive clusters from in average clusters. 
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The flux ratios v.s. the cluster mass are shown in Figure 6. Note that the SZ flux are 
represented in positive values. The flux ratio decreases with mass, since the SZ flux increases 
with M vir faster than the radio flux does. The flux ratio also decreases with redshift under 
the presently adopted no-evolution hypothesis. Moreover, at v ~ 30 GHz, the radio flux 
are higher and the SZ flux are lower than those at 90 GHz, and therefore the flux ratios of 
the two frequencies differ by a large margin. In SZ surveys for cluster mass ranging from 
M~2x 1O 14 M to 2 x 1O 15 M , the radio contamination at 30 GHz should be a serious 
concern for small clusters even at z—1 (~ 10%), whereas the contamination for small clusters 
at 90 GHz is lower than 10% even at ,2=0.25. 



5. Conclusion 

We study the CRLFs from 1.4 GHz to 180 GHz, based on the flux-limited 1.4 GHz 
observation in conjunction with studies of radio spectral index at higher frequencies. The 
lack of comprehensive radio spectra of galaxy clusters remains a problem for predicting the 
radio confusion on the cluster SZ effect. In our analysis, the cluster radio luminosity functions 
at high frequencies are obtained by converting the observed low-frequency ones with the 
observed spectral index distributions of radio sources. We obtain the frequency dependence 
of the CRLF from two spectral index studies. Assuming that CRLFs do not evolve, the 
redshift dependence of CRLFs can be derived as well. Based on this information, we calculate 
the total flux of radio point sources in a cluster for mass ranging from M — 1 x 10 14 M Q to 
5 x 1O 15 M at high frequencies and various redshifts. Apart from few exceptional clusters, 
our predictions agree with L095a (at 1.4 GHz) and CGHJC (at 28.5 GHz) observations in 
a self-consistent manner. 

We also give estimates of the confusion from radio point sources to the SZ flux. As 
demonstrated in Figure 6, the radio confusion to lower frequency {y ps 30 GHz) SZ mea- 
surements poses a severe problem for cluster of M < 1O 15 M , whereas at higher fre- 
quency {v 90 GHz) the radio confusion causes less a problem even for small clusters 
(M ~ 1O 14 M ). For observed frequencies substantially lower than 30 GHz, such as 15 GHz 
proposed to be operated at AMI, one may use Eq.(15) to estimate the expected level of radio 
confusion. Take (a) — — 1; the 15 GHz radio flux is twice higher but the SZ flux is 4 times 
weaker (due to the v 2 Rayleigh- Jeans tail) than those in the 30 GHz observation, and the 
confusion flux ratio becomes 8 times higher than that of 30 GHz. 

The estimate for the expected total confusion to a cluster of given virial mass can be 
translated to the expected loss of SZ signals, which results in reducing the cluster count 
and in turn affecting the inferred values of cosmological parameters. As has already been 
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demonstrated, the amplitude of density fluctuations averaged over 8/i _1 Mpc, <7 8 , is most 
sensitive to the SZ cluster count (Fan & Chiueh 2001). The radio confusion thus yields a 
lower value of the inferred og. Nevertheless, with the aid of a sensitive telescope for pointed 
observation, the radio confusion can partially be removed by subtraction, when the radio 
sources are sufficiently strong. However weak radio sources remain to be a problem as they 
are bound to blend into the SZ signals and cannot be subtracted away. The spectral-index 
distribution n(a) can in principle reveal the relative population of strong and weak sources 
at high frequencies. But this problem can be assessed in a different way. We have compared 
the observed radio flux of CGHJC (28.5 GHz) and L095a (1.4 GHz) with the predicted total 
radio flux in Figure 4 and 5, and the data show large scatters around the expected flux. It 
suggests that many sources are weak and can be difficult to get removed by subtraction. 
Therefore, a good SZ observing strategy is probably to set the SZ flux limit substantially 
higher than the expected flux of radio sources in a cluster. For 90 GHz observations, a mJy 
flux limit is adequate for this purpose. 

This work was supported in part by NSC 90-21 12-M-002-026 from the National Science 
Council, R.O.C. and by the National Science Foundation of China. 
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a 

Fig. 1. — Normalized spectral index distribution of radio sources. The circles are from 
SRA96 (254 sources) and the solid line is the best fit; while the crosses are from CGHJC 
(53 sources) and the dash line is the best fit. CGHJC sample is biased toward very massive 
clusters (> 10 15 M Q ) and the spectral index distribution is more irregular at low and high a 
than that of SRA96. 




Fig. 2. — Predicted cluster radio luminosity functions at v— 
lines), and at 90(l+z) GHz (the dotted lines) for 2=0.25. 
from L096. 



28.5(l+z) GHz (the dash-dotted 
The data points are at 1.4 GHz 
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Fig. 3. — Coefficients of the luminosity-weighted CRLFs [log from (a)SRA96's n(a), 
and (b)CGHJC's n(a) in frequency. 
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v Q =28.5 GHz v Q =90 GHz 




virial mass (solar mass) virial mass (solar mass) 



Fig. 4. — Total radio flux at z/=28.5 GHz, 90 GHz in cluster virial mass with redshifts 
2=0.25, 0.5, 1.0, respectively. The solid lines are converted from SRA96's n(a), and the 
dash lines are from CGHJC. Triangles represent the radio flux of 25 clusters below ,2=0.3 in 
CGHJC. 
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Fig. 5. — Total radio flux at 1.4 GHz from our predicted CRLFs. The data points are 
boosted up to £=0.09 from L095a. 
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Fig. 6. — Flux ratios at z/=28.5 GHz, 90 GHz in cluster virial mass with redshift z=0.25, 
0.5, 1.0, respectively. The curves are the same as the former. 
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Table 1: Fit parameters for the dependence of the luminosity- weighted CRLFs [log P u f u ] 
from n(a) in SRA96 and CGHJC, respectively. 



Author 


i/ GHz 


z 


v GHz 


Ci 


7i 


c 2 


72 


logP ft 


SRA96 


1.4 


0.09 


1.526 


1.75 


0.85 


33.65 


-0.43 


24.8 




28.5 


0.25 


35.625 


2.11 


0.82 


29.61 


-0.36 


23.4 




28.5 


0.5 


42.75 


2.23 


0.81 


28.95 


-0.33 


23.3 




28.5 


1.0 


57.0 


2.45 


0.80 


28.10 


-0.31 


23.2 




40.0 


1.0 


80.0 


2.76 


0.78 


27.11 


-0.27 


23.1 




90.0 


0.25 


112.5 


3.11 


0.76 


26.03 


-0.24 


23.0 




90.0 


0.5 


135.0 


3.31 


0.75 


25.31 


-0.21 


22.9 




90.0 


1.0 


180.0 


3.64 


0.73 


24.30 


-0.17 


22.8 


CGHJC 


1.4 


0.09 


1.526 


1.75 


0.85 


33.63 


-0.43 


24.8 




28.5 


0.25 


35.625 


2.11 


0.82 


31.66 


-0.42 


23.8 




28.5 


0.5 


42.75 


2.21 


0.82 


31.47 


-0.42 


23.7 




28.5 


1.0 


57.0 


2.28 


0.81 


31.30 


-0.42 


23.6 




40.0 


1.0 


80.0 


2.40 


0.81 


31.02 


-0.41 


23.5 




90.0 


0.25 


112.5 


2.52 


0.80 


30.74 


-0.41 


23.4 




90.0 


0.5 


135.0 


2.53 


0.80 


30.55 


-0.40 


23.3 




90.0 


1.0 


180.0 


2.63 


0.79 


30.29 


-0.40 


23.2 



